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ihtbocootick 

^uipakent  diYtlcpsd  in  tills  Laboratory  pwits  the  determination  of 
eig*it  of  the  dynaeic  coefficients  useful  in  describing  the  force  and 
accent  reactions  on  a  sufcnerged  bo$r  Moving  in  cater,  times  coefficients 
p  ^ise  the  partial  derivatives  of  aoaent  (about  the  yaw  aide)  and  of 
«*2oe  (in  the  horisontal  plans,  and  perpendicular  to  the  longitudinal  axle) 
with  respect  to  velocity  arJ  acceleration  coaponente  in  specified  air  mo¬ 
tion*  <  So  long  es  tbs  instantaneous  angles  of  attack  vc  snail  and  scale 
effects  are  absent,  these  coefficients  have  constant  values.  1  complete 
list  of  coefficients  is  given  in  Ref.  (1),  as  are  definitions,  sign  con¬ 
ventions  emu  formulas  for  Making  the  coefficients  nondinenaional.  The 
eight  coefficients  tabulated  below  are  those  pertinent  to  lateral  trans¬ 
lation  and  rotation  about  ths  yaw  axis  for  a  body  of  revolution: 

H  •  coefficient  of  rotary  aoaeat  derivative 
r 

Nj1  virtual  sonant  of  inertia  coefficient  (angular  acceleration) 
coefficient  of  static  sonant  derivative 

N*.*  virtual  ucnent  of  inertia  coefficient  (lateral  acceleration) 

X  *  coefficient  of  rotary  force  derivative 
r 

T*,'  virtual  inertia  coefficient  (angular  acceleration) 

X  *  coefficient  of  static  fores  derivative 

Xj*  virtual  inertia  coefficient  (lateral  acceleration) 

where  the  prim  Indicates  that  the  coefficients  are  in  dimensionless 
foes. 

It  is  the  purpose  of  tha  experimental  program  undertaken  at  tM.s 

e 

laboratory  to  determine  the  numerical  values  of  the  above  quantities  for 


2 


th*  Basic  Finner  fissile  {Fig.  1).  Because  of  the  r^ulreu  differences 
in  the  experimental  methods,  however,  the  progreat  ires  divided  into  two 
parts.  This  report  deals  only  with  Fart  1,  end  is  restricted  to  the 
following  quantities*  N^*,  X^*,  XyS  end  the  linear  combinations 

Hr*  -  end  Yp*  «  ¥*,'  .  Remaining  quantities  will  be  the  subject  of 
another  report. 


mmxw 

The  procedure  used  to  obtain  the  dynamic  moment  coefficients  is  based 
on  the  measurement  of  the  amplitude  ratio  and  phase  angle  o  the  model  while 
it  is  oscillated  at  a  known  frequency  in  the  moving  stream  of  water.  For 
this  purpose  the  Angular  Dynamic  Balance  was  used,  in  conjunction  with  the 
High  Speed  later  Tunnel  at  this  Laboratory.  A  complete  discussion  and  des¬ 
cription  of  the  equipment  is  given  in  Ref.  2,  but  a  brief  summaiy  is  given 
hires 

The  balance  can  be  regarded  as  consisting  of  several  separate  function¬ 
al  components  (Fig.  2)  •  The  first  o  these  is  the  platform  -  platform 
drive  assembly.  This  unit  consists  of  a  platform  to  which  can  be  imparted 
a  small  angular  oscillation  of  known  amplitude  and  controllable  frequency. 

To  It  is  coupled  the  spindle-model  assembly  which  is  free  to  rotate  on  ball 
bearings  and  which  is  provided  with  a  mter  seel  to  prevent  leakage  of  tunnel 
watir  to  the  outside.  The  coupling  member  is  a  calibrated  torsion  spring 
of  controllable  stiffhees.  Tiro  mirror-lens  systems,  one  of  which  is  carried 
by  the  driving  platform  and  the  other  by  the  spindle,  enable  the  observer  to 
measure  the  angular  deflections  of  each  of  these  components  with  respect  to 
ground.  By  projecting  the  image  of  a  pulsing  light  source  (triggered  by 
contacts  on  the  drive  assembly?  the  angular  deflection  at  a  given  instant 
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in  t he  cycle  can  bn  matured,  Tim  contactor  body  i*  rotated  through  known 
angular  increment*  and  the  deflection*  taken,  by  plotting  these  data  and 
by  representing  the  observed  notions  in  Fourier  series,  the  amplitude  end 
phase  angles  of  the  fundamental  motion  of  both  driver  and  spindle  can  be 
determined* 

A  self-contained  lateral  force  balance,  using  bonded  strain  gages, 
serves  as  the  center  section  of  the  model  assembly*  This  unit  measures  the 
instantaneous  sitte  force  acting  on  the  model  and  relates  it  to  the  instan¬ 
taneous  position  of  the  driving  platform* 
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Up  to  this  point  the  balance  has  been  regarded  as  a  harmonic  oscillator 
which  permits  complete  analysis  of  the  motions  of  both  the  driving  platform 
and  of  the  spindle-model  assembly*  By  running  the  balance  in  air,  and  by 
performing  the  necessary  computations,  the  spindle-model  assembly  is  found 
to  possess  a  certain  moment  of  inertia  and  little  or  no  damping  with  respect 
to  ground*  Similarly  the  internal  (lateral  force)  balance  shows  that  the 
center  of  gravity  of  the  model  assembly  coincide©  very  closely  with  the  axis 
of  angular  oscillation*  This  fact  is  indicated  by  the  absence  of  aiy 
measured  periodic  lateral  force  in  the  presence  of  the  angular  motion. 

*ben  the  tunnel  is  filled  with  water  which  is  moving  relative  to  the 
model,  the  behavior  of  the  model  is  strikingly  different.  Analysis  of  the 
motions  of  the  platform  and  spindle  shows  that  the  spindle  "sees*  a  larger 
maaa,  coupled  to  ground  through  a  spring,  and  possessing  a  degree  of  damping 
with  respect  to  ground.  This  spring  can  either  aid  or  oppose  the  model  dis¬ 
placement  and  the  damping  can  be  positive  or  negative*  In  similar  fashion, 
the  internal  balance  shows  that  periodic  sids  forces  are  present  which  behave 


u 


am  though  caused  by  oscillating  a  larger  mass#  coupled  to  ground  by  a 
spring,  and  also  possessing  sons  amount  of  deaping. 

The  relationship  between  the  instantaneous  forces  and  moments  exerted 
by  the  spindle  on  the  model  and  the  resulting  motion  of  the  model  can  be 
given  in  terms  of  the  effective  physical  parameters  of  the  system  and  the 


angular  displacement  p0 

of  the  models 
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(2) 

idlers  the  sign  convention  has  been  chosen  to  facilitate  final  expression  of 
the  desired  hydrodynamic  coefficients.  In  these  equations 

a  *  total  apparent  mass  of  the  model 

b  =  lateral  force  damping  coefficient  of  model  with  respect  to  ground 
k  »  effective  lateral  spring  rate  with  respect  to  ground 
I  a  total  effective  moment  of  inertia  of  model  in  water 
B  a  torsional  damping  moment  coefficient  of  model 
K  *  effective  torsional  spring  rate  with  respect  to  ground 

These  apparent  physical  properties  of  the  model  and  model-spindle  assembly 
are  a  result  of  the  hydrodynamic  forces  acting  on  the  model.  For  the  small 
oscillations  executed  by  the  model  and  for  the  small  instantaneous  angles  of 
attack,  the  hydrodynamic  react  -  one  can  be  regarded  as  comprising  the  sum  of 
the  indiv*  dual  reactions  arising  from  displacement,  velocity,  and  acceleration 
in  the  direction  of  motion.  Letting  N  equal  the  total  hydrodynamic  moment 
acting  on  the  body  about  the  yaw  axis  and  T  equal  the  total  hydrodynamic 
lateral  force,  then 

K  *  Hu  +  H.u  +  N  v  *  H.*  +  Hr  +  H.r 
u  u  v  v  r  r 


(3) 
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and 

X  *  Xu  +  T*u  +  If  +  T-v  +  Ir  +  X*r 
u  u  v  v  r  r 

where  subscripts  indicate  partial  derivatives  in  accordance  with  the  list 

# 

of  symbols  on  page  13  end 

U  **  u  velocity  of  bocfcr  axis  relative  to  fluid  at  infinity. 

m 

r  «  Pq  angular  velocity  of  laodal  about  the  yaw  axis. 

v  <t£-U  P  component  of  velocity  of  bodty’  axis  relative  to  fluid 
in  the'  direction  perpendicular  to  the  longitudinal 
axis  and  in  the  horixontal  plane. 

v^-U  p  component  of  acceleration  of  body  relative  to  fluid 
in  the  direction  perpendicular  to  the  longitudinal 
axis  and  in  the  horiaontal  plane. 

Since  the  body  in  question  possesses  rotation  symmetry,  the  following  terms 
must  vanish  for  constant  u. 

Nu,  H^u,  Tuu,  and  X^u  . 

Equation  (3)  can  therefore  be  rewritten  in  &* 
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For  tho  present  case,  the  center  of  gravity  of  the  test  body  was  very  nearly 
at  the  center  of  rotation,  so  that  all  side  force  reactions  were  hydrodynamic 

in  (vpiirln.  cap 

mpiS*  vw»  vest 

I  -  -  \  (5) 

JU*o*  the  gumat  reaction  on  tha  apindl.  due  to  duping  and  dtaplacernffit  in 
air  was  saroj  therefore  th.  only  nontydrocc-nuic  ruction  w  du.  to  th. 
aoatnt  of  inartt*  of  the  apindlo-ncdal  asaeably,  or 
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»  *  h  Po  "  H«  i6) 

mo  U»t  by  combining  ten*  in  (1),  (2)  end  (U)»  and  substituting  in  (5) 
and  (6),  the  hydrodynesdc  coefficients  w*n  be  identified  with  the  effective 
physical  quantities  characteristic  of  vibrating  ness  system  t 
>  *  -I.  b  *  (Tr  -  0)  k  »  -Ty  U 

(7) 

I,  =  -N-  B  ■  (N  -  H  U)  K  -  HU 

IT  T  V  f 

where  U  Is  the  constant  relative  fluid  velocity  at  a  great  distance  from 
the  bodjjr  • 

The  experimental  problem  then  reduces  to  one  of  finding  the  numerical 
valuta  of  these  effective  physical  constants,  m,  k,  B,  and  K,  by* 
whatever  means  available.  The  (frnamic  balance  and  its  accessory  instru¬ 
mentation  provides  this  means*  If  the  equations  of  motion  for  the  system 
are  formulated,  the  following  pair  of  equations  is  obtained! 


z»h  *  hK  =  Xd(?2  -  Pl>  ’  *,  *1 

I  p0  +  B  K  +  M,  a  I,  (px  -  p0) 
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» 
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where  the  symbols  are  as  defined  on  page  12  and  as  shown  in  the  equivalent 
linear  analogy  of  the  moment  measuring  apparatus  in  fig.  2.  From  these  ex¬ 
pressions  can  be  obtained  the  steafy- state  solution 
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Tv* 


Kd+Ka 


s 


I  «  icoB, 
<JS -  +  _ £ 

K  K 


cos  (  d  -  6  )  -  i  sin  (  $  -  } 


i 


t 


i 

t 

t  f 


(9) 


*n<r 


which  reduces  to 


Zt  +  t  I®2  +  I»b"| Kd+  K# 


cos  («2  -  atj_)  -  1  sin  («2  -  3^) 


J& 

Vi 


when  ten*  containing  the  negligible  quantities  I  and  B.  are  dropped. 

s  a 

The  real  part  of  this  equation  identifies  the  in-phase  component  and 
&he  imaginary  part  identifies  the  quadrature  conponent.  Fron  this  ex¬ 
pression  the  effective  physical  constants  I,  B,  and  K,  can  be  found, 
SLcoording  to  the  method  to  be  described  later.  Similar^,  the  expression 


—  cos  ($Q  -  -  i  sin  (  6q  -  d2)  m  « 

^  L  K. 

h+%a  V**  *®»  ’ 

-  -! - + - -  0C(flL-«  ) 

K.  *.  K.  T.  2 

«  B  8 


cos  -  «2)  -  i  sin  (  ^  -  «2) 


Mrsfc  reduces  to 


(10) 


K.  K .  +  X 

cos  ($0  -  $2)  -  i  sin  ($Q  -  ^2)  »  —  -f  cos(d^*d2)  -  sin  (<L  ~dL 

S  8 


and  can  be  solved  to  gtira  the  aagplitudi  and  phase  angle  of  the  model  in  terms 


lateral  force  in  related  to  the  effective  phyeic&l  constants,  a,  bs  end 
k,  through  the  stea<fcr-atate  relation 


F  jjcoe  d,  +  1  gin  d*  ]  * 


]*-[<*- 


mar)  +  1(0  b 


(U) 


where  F  la  the  peak  value  of  the  lateral  force  and  df  la  the  phase  angle  of 
this  force  relative  to  the  notion  of  the  model*  ^  < 


REDUCTION  AND  ANALYSIS  Of  DATA 


Equation  9  is  seen  to  be  a  cceplex  expression  of  the  fora 
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a  +  lb 


*  c  ♦  id 


(12) 


which  can  be  rationalised  to  yield 

c  d 

a  =  y— g  and  b*«* 

e  +  d* 


"1 - 2 

c  *  d* 


(13) 


where 
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plotting  a  against  <0  (Fig.  3)  and  measuring  the  slope  $Ag)  and 
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y-iaterespt  (1  *  E/S0)  two  of  the  equivalent  physical  constants  are  found 
in  terns  of  the  known  spindle  spring  rate.  Similarly  if  b  is  plotted 
against  m  (fig*  k)  the  slops  of  tbs  resulting  lines  is  equal  to 
An  analagous  treatment  of  tbs  lateral  faros  equation  (11)  in  Figs*  5  and  6 
yields  the  quantities  s,  br  and  kf  these  equivalent  physical  parameters 
are  then  equated  to  the  appropriate  hydrodjrnanio  coefficients  according  to 
Iq*  (7)  and  nondiaeneionalised,  as  indicated  in  Ref*  X*  the  resulting 
dynamic  coefficients  have  been  plotted  against  Reynolds  maker  in  Figs*  7 
to  12.  The  use  of  Reynolds  number  as  the  similitude  parameter  does  not 
imply  an  a  priori  judgment  that  viscosity  effects  are  paramount  in  deter* 
mining  the  magnitude  of  these  coefficients*  Indeed,  for  the  experimental 
setup  available  for  this  study,  where  the  model  size  and  fluid  viscosity 
were  fixed,  the  Reynolds  number  change  indicates  only  a  change  in  velocity. 
The  effect  of  the  presence  of  the  model  support  spindle  and  spindle 
shield  was  estimated  by  the  image  method*  That  is,  supplementary  runs  were 
mads  with  a  dussy  support  spindle  shield  on  the  side  of  the  model  opposite 
the  real  support  spindle  and  shield*  This  dummy  did  not  touch  the  model* 
The  data  points  taken  with  this  imago  in  place  are  indicated  by  solid  sym¬ 
bols  in  the  final  curves*  It  is  assumed  that  the  presence  of  the  dummy 
shield  produces  the  same  effect  as  the  original  functional  spindle  and 
spindle  shieldj  hence  the  final  curve  is  taken  to  be  displaced  an  amount 
equal  to  the  difference  between  the  original  runs  and  the  image  runs,  and 
in  the  direction  of  the  original  runs*  In  mmm  eases  the  presence  of  the 
image  strut  had  no  effect  on  the  values  of  the  coefficients!  for  these  no 
image  rune  wore  plotted* 
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shook;  or  results 

(1)  Virtual  Moment  of  Inertia  Coefficient  (Fig.  7)* 

This  coefficient  is  nearly  independent  of  velocity,  m  predicted 
by  perfect  fluid  theory,  and  shears  little  or  no  sensitivity  to  the  presence 
of  the  support  spindle. 

(2)  Mvl,  Coafficiant  of  Static  Ifcaant  Derivative  (Fig.  8)t 

This  coefficient  is  a  Measure  of  the  moment  component  due  to 
velocity  perpendicular  to  the  longitudinal  bo<fcr  axle.  For  seall  angles 
of  yaw  this  is  seen  to  be  Identical  to  the  static  moment  rats,  d  C^/da, 
and  for  this  teas  on  it  serves  as  a  good  check  on  the  experimental  and  com¬ 
putational  procedures*  Ths  agreement  between  the  values  reported  In  Fig.  8 
end  those  reported  in  Kef.  (3)  are  very  good,  even  as  to  the  effect  of  the 
presence  of  the  support  spindls* 

« 

(3)  Hr*  -  N^'.  Combined  Moment  Coefficient  (Fig*  9)* 

Actually  a  linear  coabination  of  the  coefficient  of  Rotary  Moment 
Derivative  and  the  Virtual  Moment  of  Inertia  Coefficient  (to  lateral 
acceleration),  cannot  be  resolved  without  further  experiments,  which  ere 
to  be  conducted  at  a  later  time* 

(U)  Y* »,  Virtual  Inertia  Coefficient  (angular  acceleration)  (Fig.  10)* 
r 

This  quantity  is  negligibly  small  in  the  presence  of  ths  ponderable 
value  of  Yt*  . 

(5)  IYS  Coefficient  of  Static  Force  Derivative  (Fig.  11)* 

The  analog  to  item  (2)  above,  can  )t  identified  with  the  slope  of 
the  static  lift  coefficient  curve  a  C^/da  found  by  conventional  static 
balance  measurements  *  As  with  item  (2),  the  agreement  with  previous  static 


measurements  (Ref*  3)  is  excellent* 

(6)  T  1  -  I*1,  Combined  Lateral  force  Coefficient  (Fig*  12) 
r  v 

As  with  item  (3)#  this  coefficient  cannot  be  completely 
analysed  until  the  second  phase  of  this  investigation  is  completed* 


rtw  cm  qnnttr.fi 

The  motion  of  the  tost  body  is  restricted  to  tm  plane  of  yaw,  The  sym¬ 
bols  used  to  do  scribe  this  motion  end  the  associated  hydro<frnaaic  reactions  on 


the  body  are  in  the  greater  part  Identical  to  those  rece&mended  in  the  Tech¬ 
nical  and  Research  Bulletin  V,  1-5  of  the  Society  of  Naval  Architects  and 
Karine  Engineers  titled  "Nonenclature  for  Treating  the  notion  of  a  Submerged 
Body  through  a  Fluid!1  *  These  symbols  art  marked  with  an  aatariskO)  in  the 
following  list. 

The  coefficients  have  been  made  ncn-dimensional  in  accordance  with  p*  > 


cedar es  outlined  in  the  d>ove  publication.  Dimensionless  quantities  are 
prlmsdf  all  other  quantities  are  in  the  pound-foot-sscond  system  of  units. 


*  A  cross  section  area  of  model 

Art,  A.*  peak  half  amplitude  of  notion  of  body,  support  spindle 
°  A  4  and  drive  platform 

b  damping  coefficient  in  gq,  (1) 

B  torsional  damping  coefficient  in  Sq«  (2) 

B  torsional  damping  rate  of  spindle  seal 

d  diameter  of  model 

?  peak  magnitude  of  lateral  force  in  the  horio  ntal  plane 
i  =  */T 

I  =»  +  If  total  effactirs  ament  of  inertia  of  model  in  water  (Eq.  (2)) 

1^  moment  of  inertia  of  model 

apparent  moment  of  inertia  of  model  assembly  due 
x  to  hydro<fynamic  reactione 

2m  moment  of  inertia  of  model-support  spindle 

effective  spring  rate  due  to  fcydro<fcrna»ic  reactions  (Eq,  (1)) 
effective  torsional  spring  rate  due  to  hydrodynamic  reactions  (Eq*2) 

torsional  spring  rate  of  drive  spring  which  couples  the 
driving  platform  to  the  spindle-model  assembly 

torsional  spring  rate  of  model  support  spindle 

length  of  body  - 

apparent  mass  of  model  due  to  lydrcHJtynamic  reactions  (Eq,  (1)) 


k 

X 


"d 


X. 


* 
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US*  or  SXKBOUi  (eont'd) 


*  y  yydrodynanie  moment  acting  on  the  model 

N  »  -  K  Meant  on  body  supplied  by  support  spindle 

9 

*  =  H  1  {Xjtpkf^ V)  coefficient  of  rotary  ament  derivative 

*  IL  *  M- 1  (1/2 pAd^)  virtual  moment  of  inertia  coefficient 

"  *  (ansular  ieoiluntionl^ 

*  Hv  =*  Ht*  (l/2pAd  U)  coefficient  of  static  ament  derivative 

*  N*  3  N-*  (l/2pid2)  virtual  moment  of  inertia  coefficient 

T  T  (lateral  acceleration) 

r  *  £  angular  velocity  about  the  yaw  axis 

*  0  velocity  of  txxfer  axis  relative  to  fluid  at  infinity 

*  u  component  of  velocity  of  body  axis  rslatlve  to  fluid 

in  the  direction  yewuii^i  to  the  longitudinal  eyf# 

*  v  component  of  velocity  of  body  axis  relative  to  fluid 

in  the  direction  perpendicular  to  the  longitudinal 
axis  and  in  the  borisontal  plane 

*  I  hydrodynamic  lateral  forot  acting  on  the  model 

I  9  -  7  lateral  force  on  the  aodel  supplied  by  the  support  spindle 

'S 

*  Yr  a  Ir*  (l/2pAdU)  coefficient  of  rotary  force  derivative 

*  =  T^f  (l/2pld^)  virtual  inertia  coefficient  (angular  acceleration) 

*  Tv  =»  Xf*  (l/2pAU)  coefficient  of  static  force  derivative 

*  T  •  ~  T. 1  (l/2pAd)  virtual  inertia  coefficient  (lateral  acceleration) 

V  T 

8,  pr,  angular  displacement  about  the  ym  axis  of  the  model, 

°  A  spindle,  and  drive  platform 

a  ,  d?  phase  constant  of  the  motion  of  the  model,  spindle,  and 

*  *  drive  platform 

phaate  angle  of  the  lateral  force  relative  to  model  motion 
p  mass  density  of  fluid 

«o  w  2rrf  circular  froqueacyj  f  *  frequency  in  cycles  per  second 
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Fig.  1  -  Important  geometrical  dimensions  of  the  Basic  Finner  missile. 


driving  platform  support  spindle  mooel 


K,  TORSIONAL  SPRING  RATE  OF  SPiNDwE 
B  s  TORSIONAL  DAMPING  RATE  OF  SPiNOLE  SEAL 
I  s  MOMENT  OF  INERTIA  OF  SUPPORT  SPiNDlE 
I  t  MOMENT  OF  INERTIA  OF  MODEL 
I  f  APPARENT  MOMENT  OF  1NERT1A 
*  EFFECT] /E  TORSIONAL  SPRING  RATE  DUE  10 
HOROOfNAMiC  REACTIONS 
B  EFFECTIVE  TORSION  DAMPING  RATE  DUE  TO 
HYDRODYNAMIC  REACTIONS. 
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Fig.  2 


Schematic  diagram  of  dynamic  bila  ncc  driving  platform  ami 
moment  measuring  components. 


Fig.  3  -  Moment  data  reduction  plot:  In-phase  component 
(real  part  Eq.  9). 


Fig  4  -  Moment  data  reduction  plot;  Quadrature  component 
(imaginary  part  Eq  9) 


COMBINED  FORCE  COEFFICIENT  COEFFICIENT  OF  STATIC  FORCE  DERIVATIVE 


1 


Fig.  11  -  Coefficient  of  static  force  derivative  as  a  function 
of  Reynolds  number. 


Fig.  12  -  Combined  lateral  force  coefficient  as  a  function 
of  Reynolds  number. 
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